Radial Band Structure of Electrons in Liquid Metals 
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The electronic band structure of a liquid metal was investigated by measuring precisely the evolution of angle- 
resolved photoelectron spectra during the melting of a Pb monolayer on a Si(lll) surface. We found that the 
liquid monolayer exhibits a free-electron-like band and it undergoes a coherent radial scattering, imposed by 
the radial correlation of constituent atoms, to form a characteristic secondary hole band. This unique double 
radial bands and their gradual evolution during melting can be quantitatively reproduced, including detailed 
spectral intensity profiles, with our radial scattering model based on a theoretical prediction of 1962. Our result 
establishes the radial band structure as a key concept for describing the nature of electrons in strongly disordered 
states of matter. 

PACS numbers: 71.22.+i, 71.23.-k, 73.20.-r, 79.60.-i 



Central to understanding various electronic properties of 
crystalline materials is the band theory based on periodic elec- 
tron waves. In reality, however, a huge variety of materials 
are subject to significant disorders to disturb periodic crys- 
talline structures [T-^, making it challenging to describe their 
electronic structures |5-8|. Most notoriously, the nature of 
electrons in a liquid metal (LM) without any long-range or- 
der ||3|J6l has been extremely difficult to understand and no 
well-established idea exists so far to connect their atomic and 
electronic structures. 

In recent years, important progresses were made to under- 
stand the atomic structure of LMs, establishing the local con- 
figurations of liquid fragments ll9l- [T2l and their global radial 
coiTelation 19 ^161 . In sharp contrast, while various theoret- 
ical ideas were proposed since 1960's to describe the elec- 
tronic structure of LMs lfT7l42T1 . critical experimental infor- 
mation to develop these ideas has been lacking. To our knowl- 
edge, there is only one experimental report to directly map 
out the band structure of LMs with angle-resolved photoemis- 
sion (ARPES) f6l. It used a liquid Pb monolayer supported 
on a crystalline Cu(lll) surface, where the momentum {k) 
resolved information can be preserved in the photoemission 
process |6|. This work found two persistent Fermi surfaces 
(FSs) to have very short electron coherence lengths in LMs. 
However, the origin of these two bands and their markedly 
distinct coherence lengths have been mysterious. In the ex- 
perimental point of view, this limitation was partly due to (/) 
a too high melting temperature (T^ = 568 K) to resolve the 
fine FS shape and (//) the intermixed signal from the metallic 
substrate. 

In this Letter, we overcome such limitations by choosing a 
Pb monolayer on a semiconducting Si(lll) substrate with a 
band gap, which has anomalously low T^. ~ 270 K ll22l - E4l 
owing to its intrinsic diffusion property 1251 l26l . This low 
allows us to trace the fine evolution of ARPES spectra in the 
melting process at a relatively low temperature with reduced 
thermal broadening. We identify the electronic fingerprint of 
coherent radial scattering, which is a missing link to relate the 
atomic and electronic structures in LMs. Our radial scattering 



model based on an old theoretic prediction ifTTI explains the 
origin of unique double radial bands as well as their mysteri- 
ous coherence lengths, providing comprehensive understand- 
ing of the electronic band structure of LMs. 

The n-type Si(lll) wafer was thermally cleaned in the 
ultra-high vacuum chamber (6.5 x lO^'^^ torr). The Pb mono- 
layer was fabricated by in-situ thermal deposition of Pb onto 
the Si( 11 1)7x7 surface held at 500 K. The melting transition 
is completely reversible and does not depend on the annealing 
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FIG. 1: (color online), (a) Crystalline structure of the Pb mono- 
layer on Si(lll) with the areal density of 9.4 Pb atoms nm^^ |28|. 
Dashed lines represent the \/7 x \/3 surface unit cell, (b) LEED pat- 
terns taken at marked temperatures with an electron energy of 8 1 eV. 
Dashed lines indicate Si(l 1 1)1 x 1 spots, (c) Intensity and correlation 
length as a function of temperature around Tc obtained by the peak- 
fitting analysis of a \/7 x \/3 superspot in (b) using the correlation 
function oc exp(-|r|/^), r — !> oo 1 15 1. The solid line is a fit with the 
critical exponent of 0.37 predicted in a 2D liquid theory 13J. 
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FIG. 2: (color online). Experimental FS data at various marked 
temperatures. Overlaid in 90 K data is the band theory calculaton 
1271 1281 and dashed lines indicate the Si(ll 1)1x1 Brillouin zone. 
Across Tc ~ 285 K, the complex FS reduces into a large outer and 
a smaller inner FSs (white lines), the latter of which gradually trans- 
forms from circle to hexagon as further increasing temperature. The 
other weak features (grey lines) are due to remaining substrate poten- 
tials. The white dots at 600 K are band-peak positions obtained by 
fitting the radial spectral profiles. The red dots show the midpoints 
of inner and outer FSs in the radial direction. 
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FIG. 3: (color online), (a) Band dispersions taken along the arrow in 
Fig. 2. (b) Typical k distribution curve taken from (a) at the Fermi 
energy. The red line overlaid is the result of the model calculation in 
Fig. 4. Temperature dependence of (c) electron coherence length and 
(d) Fermi wavevector of the outer (closed circles) and the inner (open 
circles) bands. The coherence length is obtained from the spectral k 
width. The Fermi wavevectors of inner and outer bands are offset 
respectively by ±0.325 A^^. Red dots in (d) show the center of 
inner and outer Fermi wavevectors (radial zone boundary). 



history up to 700 K. The density of metal deposits was thor- 
oughly calibrated with the well-known phase diagram |23|. 
ARPES spectra were taken using a Scienta R4000 electron 
analyzer and a high-flux He discharge lamp (He la, 21.2 eV). 
The energy and angular resolutions were 10 meV and 0.25°, 
respectively. 

Figure 1(a) shows the crystalline structure below T^- The 
Pb layer has two groups of atoms, one sitting on top of surface 
Si atoms (blue, Tl or Tl') and the other at the hollow sites in 
between (red, H3) |l23]|3j. These result in the ^7 x ^3 peri- 
odicity [dashed lines in Fig. 1(a)] with respect to Si(l 11)1x1. 
The melting of the Pb monolayer is demonstrated in the tem- 
perature series of diffraction patterns [Fig. 1(b)]. In clear con- 
trast with the Si-related 1x1 spots (dashed lines), the higher 
order superspots from the Pb monolayer diminish rapidly at 
around room temperature. The quantitatively measured inten- 
sity of Pb-related spots in Fig. 1(c) reveals more details of the 
temperature dependence. The correlation length (^) extracted 
from the diffraction width [15] shows an abrupt downturn to 
reach to a few A. This value is on the order of the interatomic 
spacing, which manifests the liquid state. The curve fitting of 
^(r) with the critical exponent value of 0.37 predicted in the 
2D melting theory [SJ follows well the experimental results 
[solid line in Fig. 1(c)], and determines the of 285.5 K. 

Figure 2 shows the temperature series of ARPES intensity 
maps at the Fermi energy. The intensity map taken at 90 K, far 
below Tc, represents the FS of the solid Pb monolayer ll27l . 
which can be well reproduced by the band theory calculation 
(overlaid lines) [28]. It consists of two metallic bands SI and 



S2 woven complexly according to the \/7 x -\/3 periodicity 
within the silicon band gap, which originate from in-plane Pb 
6p orbitals I27u28i. 

As approaching Tc at 280 K, the complex FS pattern grad- 
ually disappears and is dominated at 290 K by two simple 
circular FSs, a large outer and a small inner one (white lines). 
Very similar FSs were also observed on Pb/Cu(lll) [6], sug- 
gesting that this double-band structure is characteristic to the 
liquid Pb monolayer. Upon further heating to 300-360 K, 
these FSs undergo reversible shape changes, gradually from 
circles to rounded hexagons, which is more obvious for the 
inner one. This finding, which was overlooked in the pre- 
vious study 1 6 1, gives an important insight into the mecha- 
nism behind the double bands as discussed below. There are 
other, much weaker, features of triangles and hexagons (grey 
lines), which are due to the weak 1x1 substrate potential [Fig. 
1(b)] ||29| . These extra features decay gradually into the back- 
ground at higher temperatures, whereas the double hexago- 
nal FSs, associated with the liquid monolayer, become more 
dominant in 600 K data, indicating their distinct origins. The 
observed double-hexagon topology of FSs, along with its ap- 
parent shape change with temperature, is not compatible at all 
with any periodic band model of a crystalline solid, calling for 
a new description beyond the scope of Bloch waves. 

Figure 3(a) shows the band dispersions taken along the ar- 
row in Fig. 2, which indicate that the outer and the inner FSs 
are electron and hole pockets, respectively. The outer band 
with a dominant intensity disperses following a free-electron- 
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like parabola, and can be traced to nearly free Pb-valence elec- 
trons with limited coherence ji6ji7j, which is a reasonable ap- 
proximation for electronic states away from the band gap in 
the theoretical model of LMs |fT9]j2TI . Supporting this idea, 
the measured electron density of the outer band (« = 2.98 x 
10^^ cm~^) amounts to the whole free-electron density of the 
Pb monolayer This, however, leaves the origin of the inner 
band mysterious. One important finding is that the outer and 
the inner bands disperse symmetrically; these two bands can 
exactly be overlapped by folding at the middle k point [dashed 
line in Fig. 3(a)]. This denies the previous assumption that 
these two bands originate from different orbitals ||6|. The band 
folding is typical for electron waves modulated by a given pe- 
riodic or aperiodic structural order |7i|8J, and suggests some 
structural order of LMs. 

Motivated by a theoretical idea ifTTl , we focus on the re- 
maining structural order in LMs, the radial atomic correlation, 
which has been well-established with diffractions ll9l [l3HT5l . 
This early theory considers simple plane waves scattered by 
screened ionic potentials. The range of ionic order {R) versus 
electron mean-free path (A) determines the qualitative nature 
of electron scattering. There are two extreme cases: When 
R is much larger than A, the situation is similar to polycrys- 
tals or glassy metals |17 |, in which many tiny crystalline do- 
mains are distributed in random fashion. The electrons would 
be globally incoherent yet locally close to Bloch states. In 
this case, the primary band would be translated according to 
the local crystalline periodicity but the lost orientational order 
distributes these translation vectors uniformly in all directions 
as in Fig. 4(a). These uniform translations of the primary 
FS would overlap to yield the secondary FS-like feature or a 
broad hole pocket. The other extreme case is an ideal liquid, 
where R is comparable (or shorter) to A ifTTl . In this case, 
the electrons are subject to a radial potential due to the radial 
atomic correlation, leading to scattering preferentially along 
the radial direction. Thus, the primary band would be radi- 
ally scattered to create a more clear hole band as in Fig. 4(b), 
which constitutes the characteristic double-band structure. Al- 
though the distinction between these two cases in reality may 
not be too obvious, the fundamental difference is the electron 
scattering nature, uniform or radial scattering. For a LM at 
a finite temperature, the electron state would be in between 
these two extremes. Moreover, as R tends to decrease with in- 
creasing temperature, the band evolution above Tc would be 
a crossover in between these two extremes with the gradually 
enhanced radial scattering. 

Based on the above idea, we performed a simple model cal- 
culation II29I . In our model, the enhancement of the radial 
scattering is parameterized by the power (p) value that deter- 
mines the angular decay of scattering intensity from the radial 
direction [Fig. 4(c)]; uniform scattering at /:> = and strong 
radial scattering at higher p values. We determined the pri- 
mary band shape from the curve-fitting analysis of the outer 
FS at 600 K (white dots in Fig. 2). This model success- 
fully reproduces the major experimental observations [Fig. 
4(d)]; the characteristic double-band structure, and more im- 




FIG. 4: (color online). Schematic illustrations of the electron scatter- 
ing and expected FSs for two models in Ref. 1 17 1, (a) the polycrys- 
talline model and (b) the ideal liquid model. The arrows represent 
examples of possible scattering vectors and the open circles are ex- 
perimental band positions at 600 K. (c) Angular distribution of scat- 
tering intensity (in log scale) from the radial direction used for (d). 
The degree of radial scattering is parameterized by the p value, (d) 
Calculated FSs based on the simple model considering the gradually 
enhanced radial scattering [that is, varying the p value in (c)] . 



portantly the gradual shape change of the inner FS from circle 
at p = 0.6 to hexagon atp = 2.0. This remarkable agreement 
between the model and the experiment corroborates the valid- 
ity of our interpretation and the old theory 1 17 1; the hole band 
comes from the radial scattering of the primary band. 

The gradual shape change of the inner FS is a consequence 
of enhanced radial scattering, indicating the continuous elec- 
tronic evolution above from incoherent Bloch waves to ra- 
dial waves. This can be related to continuous melting that 
passes through an intermediate glassy state lil4l IT6l . The 
hexagonal shape of the primary FS at a higher temperature, 
which is not very clear at the intermediate state due to over- 
lapped 1x1 bands (Fig. 2), was also consistently found in 
liquid Pb/Cu( 111) 1 6 1 . This reflects a loc al point symmetry of 
liquid Pb fragments, another important structural characteris- 
tics of LMs li9 M2l . While the hexagonal symmetry itself can 
be intrinsic to the liquid Pb monolayer |9 12 1, their orienta- 
tional alignment would be preferentially along the underlying 
substrate potential mesh [11, 13 1. This allows the observation 
of the apparent band structure from LMs 1(6, i9J. Even under 
the finite substrate influence, the electron wave only coher- 
ent over a few atoms [Fig. 3(c)] represents a strongly disor- 
dered electronic system, and would be indistinguishable from 
an ideal 2D liquid |6 |. 

The electron coherence length as a function of temperature 
is displayed in Fig. 3(c). Above Tc the inner band has a 
significantly, about 1.5 times, shorter coherence length (that 
is, broader spectral width) than the outer one. This is more 
clearly manifested in the k distribution curve in Fig. 3(b). A 
consistent observation was made for the liquid Pb/Cu(lll), 
where the authors interpreted it as a different degree of lo- 
calization and argued the distinct orbital symmetry and/or the 
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wavelength of two bands as possible origins fSl. However, 
according to our interpretation the inner band is just a pro- 
nounced trace of complex scatterings from the outer primary 
band, and its broader spectral width is natural and intrinsic to 
the scattered origin. Indeed, our model calculation does re- 
produce the distinct spectral widths and even the strange val- 
ley intensity in between [the red line in Fig. 3(b) and see 1291 
for the y direction]. In addition, the same procedure can be 
successfully extended to all binding energies throughout the 
experimental dispersions [29 1. That is, our radial scattering 
model provides a description for the whole spectral function 
ofLMs. 

For the ideal liquid model in Fig. 4(b), one can introduce 
the concept of radial Brillouin zone [17,il8J, as counterparts 
of periodic Brillouin zone in crystalline solids IJI and 'quasi- 
Brillouin zone' in quasicrystals [8 1. Connecting the midpoints 
of observed inner and outer bands in every radial direction, 
one can experimentally trace the distribution of the reference 
k points for the radial scattering. The measured midpoints are 
surprisingly isotropic (radial) as shown by red dots in Fig. 2. 
This is in clear contrast to the conventional Brillouin zone of a 
simple polygon for a 2D solid [5|. We term the double-radial 
band and the radial Brillouin zone as the radial band struc- 
ture, which is the key concept for the electronic structure of 
LMs. This would have wide implications not only in devel- 
oping more sophisticated theory but in understanding exotic 
electronic properties of various strongly disordered systems, 
for example, electron localization in glassy metals |[Tl|2l. 

From the k position of the radial Brillouin zone {\q\l2), one 
can derive the radial scattering amplitude (|^|) and the aver- 
age interatomic spacing (a = 2TTl\q\) [10|. In Fig. 3(d), we 
show the temperature dependence of radial zone boundary. It 
decreases with increasing temperature and converges to about 
1.01 A-\ which yields \q\ = 2.03 A"! and a = 3.10 A. In 
fact, this is fully consistent with the known values for liquid 
Pb/Ge(lll) in X-ray diffraction (a = 3.06 A) |[l3l and our 
own analysis for Pb/Cu(lll) (a = 3.11 A) |^|. The remark- 
able agreement on three independent systems with different 
surface lattice constants (3.61-4.00 A) points to the existence 
of universal interatomic spacing for liquid Pb fragments irre- 
spective of substrates, indicating dominant interaction within 
a 2D liquid. The radial Brillouin zone is thus a powerful con- 
cept to link atomic and electronic structures of LMs. 
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